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Abstract: We argue that by measuring higher moments of the net proton number fluctuations in heavy ion collisions 
(HIC) one can probe the QCD chiral cross over transition experimentally. We discuss the properties of 
fluctuations of the net baryon number in the vicinity of the chiral crossover transition within the Polyakov 
loop extended quark-meson model at finite temperature and baryon density. The calculation includes non- 
perturbative dynamics implemented within the functional renormalization group approach. We find a clear 
signal for the chiral crossover transition in the fluctuations of the net baryon number. We address our theo- 
retical findings to experimental data of STAR Collaboration on energy and centrality dependence of the net 
proton number fluctuations and their probability distributions in HIC. 
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1. QCD phase diagram 

Although, the important question addressed in QCD on the existence of a true 2nd order phase transition at finite 
chemical potential (critical point) has not been answered yet, nevertheless, there is recently an essential progress 
in the quantitative description of the QCD phase diagram. The lattice QCD (LQCD) provided a final value for 
the chiral cross over transition temperature at vanishing chemical potential [1]. The LQCD has also provided 
arguments that at small n ~ the chiral cross over line is the pseudo critical line of the 2nd order chiral phase 
transition belonging to the universality class of 3-dimensional, 0(4) symmetric spin models [2]. Based on the 
universality arguments the LQCD has also provided the curvature of the chiral transition line, T c (/j, c ) [2]. These 
results confirm that, at least at small values of the baryon chemical potential /i, the chiral cross over transition 
appears in the near vicinity to the chemical freezeout line [3] obtained from the analysis of particle yields measured 
in HIC. 

The numerical coincidence of thermal parameters for the chiral transition and the freezeout conditions indicates 
that the hadron resonance gas partition function, which describes chemical equilibration of particle yields in HIC, 
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Figure 1. The temperature dependence of kurtosis R.4 % 2 '■— 9X4 /X2 an( ^ ^ nc hig ner order, Xq /X2 ana - Xs I X2 > ra tios of 
cumulants for different fJ, q /T calculated in the PQM model within the FRG approach [6, 8]. The T pc is the pseudo- 
critical temperature obtained in the model at the physical pion mass. 



should describe also the QCD thermodynamics up to a near vicinity to the transition to a quark gluon plasma 
phase. Indeed, the equation of state calculated on the lattice as well as other thermodynamical observables in the 
hadronic phase, were shown to be very well quantified by the hadron resonance gas (HRG) partition function [4]. 
Already at vanishing chemical potential, i.e. under conditions realized in the high energy runs at RHIC or LHC, 
the question arises to what extent a refined analysis of freeze-out conditions can establish the existence of a chiral 
phase transition. In this lecture we will argue that even at \xb /T ~ the net baryon number fluctuations and 
their higher moments can be used to identify the chiral cross over transition experimentally [5-8]. 

2. Charge fluctuations and the chiral cross over transition 

Due to remnants of 0(4) criticality related with the chiral phase transition observed in LQCD, the free energy 
(/) near the chiral phase transition temperature T c may be represented in terms of singular (f 3 ) and regular 
contributions (f r ) as 

f(T,n q ,m q ) = f 3 (T,fi q ,m q ) + f r (T, fj, q ,m q ) , (1) 

where in addition to the temperature T we also introduced explicit dependence on the light quark chemical 
potential, fj, q = /j.b/3, and the (degenerate) light quark masses m q = m u = raj. The singular part of the free 
energy may be written as [11] 

f s (T,» q ,h) = h 1+1/S f s (z) , z = t/h 1/ps (2) 

with /3, 8 are critical exponents of the 3-dimensional, 0(4) universality class and t = ^ T y Tc + n q (^?-) 2 j, 
h s -t-tf 2 -- Here T c is the phase transition temperature in the chiral limit and to, ho are non-universal scale 
parameters. The proportionality constant K q ~ 0.06, has recently been determined from a scaling analysis in 
(2+l)-flavor QCD [2]. The scaling function f s and its derivatives have recently been calculated using high 
precision Monte Carlo simulations of the 3-dimensional 0(4) spin model [12]. 

We want to focus here on properties of moments of net baryon number fluctuations, which are obtained from 
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Figure 2. Left-hand figure: the ratio of quadratic fluctuations and mean net baryon number ((T 2 /M), cubic to quadratic (S<r) 
and quartic to quadratic (ftcr 2 ) baryon number fluctuations calculated in the HRG model on the freeze-out curve 
[5] and compared to results obtained by the STAR Collaboration [14]. The dashed curves show the approximate 
tanh^s/T 7 ) result for ncr 2 and Scr, respectively. Middle figure: The probability distributions, uncorrected for event- 
by-event counting efficiency, for the net proton number for different ccntralitics taken by STAR Collaboration in Au-Au 
collisions at ^/sjvlv" — 200 GeV [14]. The lines are the Skellam distributions calculated within HRG model [7]. Right- 
hand figure: Mean (M), variance (tx), skewness (S) and kurtosis (k) of the net proton number calculated from the 
probability distributions shown in the middle figure. 



Eq. 1 by taking derivatives with respect to f±B = [Ib/T, 



b i d"//r 4 l d"/ r /r 4 l d"f s /T 4 _ B B 

Xn 3 „ g .„ 3 „ g„ Xn,r + Xn,s ■ \?) 

In the hadronic phase and away from transition temperature the regular part Xn,r should be well described by 
the hadron resonance gas partition function which will be a reference for critical fluctuations coming from the 
singular part Xn s- Thus, any deviations from the regular i.e. hadronic gas contribution could be an indication of 
criticality due to remnants of the chiral 0(4) transition. 

Higher order moments will become increasingly sensitive to the singular part of the free energy. From Eq. 2 it is 
apparent that these moments show a strong quark mass dependence in the vicinity of the critical temperature, 



B 

Xn.s ~ < 



-{2n q ) n ' 2 h ( - 2 - a - n ^^ 5 /i n/2) (z) , for n q /T = 0, and n even 

(4) 

-(2k,)" (^) n / l (2 - Q -" )//34 /i n) ( 2 ) , for fiq/T > 



where we used 2 — a — 138(1 + 1/8). As a — —0.2131(34) is negative in the 3-dimensional, 0(4) universality class, 
the 4th order moments of the net baryon number fluctuations do not diverge yet in the chiral limit at the chiral 
transition temperature, z = 0. The first divergent moment is obtained for n = 6 if L-iq/T — and for n = 3 if 
fi 9 /T > 0. 

In the hadron resonance gas (HRG) the regular part of the fluctuations Xn,r can be directly calculated from the 
thermodynamic pressure following Eq. 3. The HRG is a mixture of ideal gases of all particles and resonances, con- 
sequently the thermodynamic pressure exhibits a factorization of T and [ib/T dependence. Under the Boltzmann 
approximation the pressure in the HRG, P HRG {T, hb) — /(T) cosh(/is/T), where f(T) contains contributions 
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from all baryons and baryonic resonances. With such p HRG there are particular properties of ratios of cumulants, 

B B B B 

X2n,r | 1 X(2n+l),r I .. A(2n + l),r, , , X2n,r , , /rr \ /r-\ 

B \HRG = 1 , g \HRG = 1 , g \HRG — tSXlh^B / i ) , — | HRG — COth(/i B / 1 ) , (5) 

X2,r Xl,r X2,r Xl,r 

which are independent of the number of baryons, their masses, degeneracy factors or decay widths. 
The above structure of cumulants ratios observed in the HRG will be modified if the singular part is included in 
Eq. 3. The Xn,a are increasingly sensitive to the order of cumulants. Thus, ratios of cumulants with different n 
should be strongly varying functions of T and /j,b when approaching a chiral cross over transition. Consequently, 
the observed deviation from that expected in Eq. 5 could be considered as a signature of the singular part 
contribution to the overall fluctuations thus, also of the chiral cross over transition. 

The generic structure of ratios of different cumulants near the chiral cross over transition is shown in Fig. 1. 
These ratios were calculated in the Polyakov loop extended quark-meson (PQM) model at the physical pion mass, 
applying the functional renormalization group (FRG) method. In FRG approach one includes quantum and 
thermal fluctuations which are needed to preserve the universal scaling behavior of physical quantities expected 
in the 0(4) universality class. In the low temperature phase the ratios X2n,r/X2,r = 1 as expected in the hadron 
resonance gas from Eq. 5. For n — 2, the kurtosis f?4 i 2 = 9x?,r/X2,r is not affected by the chiral critical dynamics 
since there is no contribution of the singular part to the first four moments as seen in Eq. 4. The observed in Fig. 
1 drop in i?4,2 is due to "statistical confinement" property of the PQM model [13]. Such behavior of kurtosis was 
first observed in LQCD calculations and was interpreted as being a signature of deconfinement in QCD [11]. 
The large deviations of (i?4,2/9) from unity in hadronic phase, which are increasing with n/T, are due to a singular 
contribution to the x? an d X2 ratio (see Eqs. 3 and 4). With increasing order of cumulants and the value of 
the chemical potential their ratios are dominated by the singular part already deeply below the chiral cross over 
transition temperature. Such behavior could be observed in HIC if freezeout appears near the chiral transition. 
Recently, the first data on charge fluctuations and higher order cumulants, identified through the net-proton 
fluctuations, were obtained by the STAR Collaboration in Au- Au collisions at several collision energies [14]. To 
explore possible signs of chiral criticality and a cross over transition, the STAR data on the first four moments 
are compared in Fig. 2 to HRG [5, 7] following Eq. 5. Different ratios of cumulants of the net proton number 
can be directly connected with measured mean (M), variance (<r), skewness (S) and kurtosis (k) [5]. 
The basic properties of measured fluctuations and ratios of cumulants are consistent with that expectated in the 
HRG model [5, 7]. This indicates, that moments of the net proton number are of thermal origin with respect to 
the grand canonical ensemble and that they freezeout along the same chemical freezout line as particle yields, close 
to the chiral cross over line. However, already such first comparison of the HRG model with STAR data reveals 
that there are deviations [5, 7, 9]. This is seen in Fig. 2 on the level of different ratios of cumulants as well as by 
comparing directly the measured probability distributions with Skellam distribution expected in the HRG for the 
net proton number [7]. The HRG model, as seen in Fig. 2, results in a broader distribution than observed in data, 
particularly at the most central collisions. This implies deviations of data on different moments from the HRG 
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model which are increasing with centrality. Shrinking of measured widths of the net proton number distribution 
relative to the HRG results, seen in Fig. 2, is to be already expected due to deconfinement properties of QCD 
[7]. We have to stress, however, that the experimental net-proton distributions in Fig. 2 are not corrected for 
event-by-event proton/anti-proton counting efficiency. While the HRG lines in this figure used efficiency corrected 
mean multiplicities from the same experiment. If uncorrected data are used in the Skellam distribution, then the 
agreement of the HRG model and measure probability distributions is found to be much better [10]. 
From the model calculations in Fig. 1, as well as from Eqs. 3 and 4, it is clear that contributions of the 0(4) 
singular part to fluctuations increase with the order of cummulants. Preliminary data of STAR Collaboration 
on Xe 1X2 ratio, taken in Au-Au collisions at the top RHIC energy, show a strong deviation of this ratio from 
the HRG model expectations with a rather moderate deviations of the lower order cumulant ratios. Recently, 
the STAR Collaboration has also observed that deviations in central Au-Au collisions of \4 lx% ratio from the 
HRG is non-monotonic in energy. Such behaviors of ratios of cumulants are to be expected due to remnants of 
criticality related with a chiral cross over transition. It is a further challenge to understand and quantify the 
observed energy and centrality dependence of these deviations from the HRG results. 
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